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Superfluid properties of the inner crust matter of neutron stars, formed by nuclear clusters im- 
mersed in a dilute neutron gas, are analysed in a self-consistent HFB approach. The calculations 
are performed with two pairing forces, fixed so as to obtain in infinite nuclear matter the pairing 
gaps provided by the Gogny force or by induced interactions. It is shown that the nuclear clusters 
can either suppress or enhance the pairing correlations inside the inner crust matter, depending on 
the density of the surrounding neutrons. The profile of the pairing field in the inner crust is rather 
similar for both pairing forces, but the values of the pairing gaps are drastically reduced for the 
force which simulates the polarisation effects in infinite neutron matter. 



I. INTRODUCTION 

The possibility of nuclear superfiudity in neutron stars 
was suggested long ago 0,0: before the first pulsars were 
actually observed. The first experimental fact pointing 
to the nuclear superfluidity in neutron stars was the large 
relaxation times which follow the sudden period changes 
( so-called "glitches'") of pulsar rotation. Later on, the 
mechanism of the glitch phenomenon itself was directly 
related to the rotational behaviour of the inner crust 
superfluidQ- Thus, according to the present models the 
glitches are generated by a catastrophic unpinning of the 
superfluid vortex lines from the nuclei immersed in the 
inner crust of neutron stars. 

Apart from the glitches phenomena, the superfluid 
properties of the inner crust matter have also important 
consequences on the cooling of neutron stars |4J. Thus, 
the heat diffusion from the interior of the star to the 
nuclear surface can be strongly suppressed in the crust 
region due to the energy gap in the excitation spectrum 
of the inner crust superfluid. 

A fully microscopic calculation of pairing properties of 
inner crust matter should be based on a bare nucleon- 
nucleon force and should account for the polarisation ef- 
fects induced by the nuclear medium. However, this is a 
difficult problem which is not yet completely solved even 
for the pure neutron matter [fj. In addition to the un- 
certainty related to the appropriate pairing force, for the 
inner crust one has also to face the problem of calculating 
the pairing properties of a non-uniform nuclear system, 
in which the density is changing significantly from the 
center of the nuclear clusters to the dilute neutron mat- 
ter in which the clusters are immersed. 

The first calculations of the inner crust matter su- 
perfluidity were done by employing semiclassical pairing 
models based on a local density approximation [|| . Later 
on the calculations were performed also in the Hartree- 
Fock-Bogolibov (HFB) approach Q. In these works the 
mean field of the inner crust matter was fixed to a Woods- 
Saxon form and the pairing field was estimated by using 
a Gogny force y| and a bare force (Argonne). The HFB 
calculations showed that the semiclassical pairing models 



overestimate the influence of the nuclear clusters on the 
pairing properies of the inner crust matter. 

In this paper we present the pairing properties of the 
inner crust matter predicted by fully self-consistent HFB 
calculations in which both the mean field and the pair- 
ing field are calculated starting from effective two-body 
forces. In this way, the effects induced upon the pairing 
correlations by the strong local variations of the parti- 
cle density and of the nucleon effective mass inside the 
inner crust matter are taken into account consistently. 
Moreover, the self-consistency allows us to study in more 
details how the properties of the nuclear clusters are af- 
fected by the external neutron gas in which they are em- 
bedded. 

In the present HFB calculations the pairing correla- 
tions are evaluated by using a density-dependent contact 
force. The calculations are done with two sets of pa- 
rameters, fixed to reproduce the superfluid properties of 
neutron matter given either by a Gogny force or by micro- 
scopic calculations which take into account polarisation 
effects [iGji- 

The article is organized as follows. In Section II we 
present shortly the structure of the inner crust matter 
used in this study and we introduce the HFB approach. 
The results of the calculations are presented in Section 
III. In the first part of this section we analyse the density 
and the mean field distributions and in the second part 
we discuss the pairing properties. The summary and the 
conclusions are given in Section IV. 

II. INNER CRUST MATTER AND THE HFB 
APPROACH 

According to standard models the inner crust matter is 
formed by a lattice of neutron-rich nuclei immersed in a 
sea of unbound neutrons and relativistic electrons ■ I n 
the density range from neutron drip density (pd w 1.4 x 
10~ 3 /9o, where /9o=0.16 fm -3 is the saturation density) to 
about half the saturation density, the nuclear clusters are 
most probably spherical and the unbound neutrons are 
in the 1 S'o superfluid phase. At higher densities, before 
the nuclear clusters are dissolved in the uniform matter 
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of the core, other non-spherical nuclear configurations 
(e.g., rods, plates, tubes, bubbles) can be formed. 

One of the first microscopic studies of the properties 
of inner crust matter was done in Ref.^2]. In this ref- 
erence the distribution of the baryonic matter in the in- 
ner crust was evaluated by using a density functional 
suggested by the density matrix expansion method and 
adjusted on the neutron equation of state of Siemens- 
Pandharipande[13j. In these calculations the pairing cor- 
relations and the spin-orbit interaction for the neutrons 
were neglected. The calculations were performed in the 
Wigner-Seitz approximation, i.e., the inner crust mat- 
ter is replaced by a set of non-interacting cells, each cell 
containing in its center a nuclear cluster surrounded by a 
gas of unbounded neutrons. In addition, in each Wigner- 
Seitz cell are uniformly distributed a number of relativis- 
tic electrons equal to the number of protons contained in 
the nuclear cluster. 

In this study we analyse the superfluid properties of the 
inner crust matter starting from a set of such Wigner- 
Seitz cells determined in Ref^l- More precisely, for a 
given total baryonic density and proton number we cal- 
culate in the HFB approach what are the mean fields and 
the pairing fields of the nuclear matter distributed in the 
Wigner-Seitz cells. 

The HFB equations are solved here in the coordi- 
nate space. The calculations are performed by using a 
Skyrme-type force for the particle-hole channel and a 
zero range pairing force in the particle-particle channel. 
In this case the HFB equations are local and for a spher- 
ically symmetric system they reduce to a set of radial 
equations [T3.IT5| : 



h(r)-X A(r) \(U i (r)\_ (U i (r) 
A(r) -h(r) + X \Vi(r) J-^y V t (r) 



(1) 



where Ui, Vi are the upper and lower components of the 
radial HFB wave functions, A is the chemical potential 
while h{r) and A(r) are the mean field hamiltonian and 
pairing field, respectively. They depend on particle den- 
sity p(r), abnormal pairing tensor n(r), kinetic energy 
density r(r) ans spin density J(r) defined by: 



Kr) = ^I>ji + l)V7(rMr) (2) 

i 
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The summations are over the whole positive-energy 
quasiparticle spectrum of the system. For the unbound 
quasiparticle states the summations should be replaced 
by integrals over the energy |l6j |. The general expres- 
sions of the mean field and pairing field on the densities 
are given in Ref . . 

The HFB equations are applied here to the inhomo- 
geneous matter of a spherical Wigner-Seitz cell, i.e., a 
nuclear cluster immersed in a sea of unbound neutrons. 
Thus, at relatively large distance from the nuclear cluster 
the particle density and the mean field should have con- 
stant non-zero values corresponding to a uniform neutron 
gas. In order to get this physical situation we impose at 
the cell radius R c the following boundary conditions for 
the HFB solutions [l^: i) even parity wave functions 
vanish at r = R c ; ii) first derivatives of odd-parity wave 
functions vanish at r — R c . With these mixed boundary 
conditions at the cell border the continuous quasiparticle 
spectrum is discretized and all the HFB wave functions 
corresponding to bound and unbound states are normal- 
ized in the Wigner-Seitz cell. 

In the HFB calculations we use for the particle-hole 
channel the Skyrme effective interaction SLy4 ^3 • This 
interaction has been adjusted to describe properly nuclei 
with a large neutron excess as well as the properties of 
neutron matter. Due to these constraints one expects 
that the nuclear matter distribution provided by SLy4 
is not far from the one estimated in Ref. Il3| . Thus, as 
seen from the calculations presented in Ref. IS] the force 
SLy4 predicts for the inner crust matter almost the same 
maximum density and proton fractions as the density 
functional employed in Ref. [T^| . 

The pairing field is calculated here with a density- 
dependent contact force of the following form 0] : 



V(r - r') = V [l - v(-) a ]5(r - r') = V ef f (p(r))S(r - r'). 
Po 

(6) 

With this force the pairing field is local and is given by: 



A(r) = V eff (p(r))K(r). 



(7) 



dr r 



Since the pairing force has a zero range the HFB calcula- 
tions should be performed with an energy cut-off for the 
quasiparticle spectrum. 

For the parameters of the pairing force we use the fol- 
lowing values: V =-430.0 MeV fm 3 , 77=0.7, and a=0.45. 
With these values and with a cut-off energy equal to 60 
MeV we reproduce approximately the pairin g pr operties 
of neutron matter given by the Gogny force |2fJ • In 
subsection 3.2 we perform also a HFB calculation with 
a second set of parameters for the pairing force adjusted 
to reproduce the gap values of neutron matter predicted 
by calculations which take into account in-medium po- 
larisation effects. 
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III. HFB PROPERTIES OF INNER CRUST 
MATTER 

As we have mentioned in the previous section, the HFB 
calculations will be performed for a set of representative 
Wigner-Seitz cells determined in Ref . [T^| . 

The maximum density for the inner crust considered 
in Ref. .12] is p ma x=0-0789 fm -3 and corresponds to a 
cell formed by 32 protons and 950 neutrons. Above this 
density the energy per baryon becomes close to the value 
of the uniform neutron system and other non-spherical 
configurations might be formed. The density region be- 
tween the neutron drip density pd and p m ax is divided 
in Ref. 01 m t° H domains. What is remarkable is that 
the cells correspondings to densities smaller than p max 
contain only 50 or 40 protons. 

From all the cells with 50 protons we will consider here 
only two representative cells, i.e., one with the density 
/9=0.0204 fm -3 and containing 1750 neutrons, and one 
with the density p—0. 00373 fm -3 and having 900 neu- 
trons. These cells have the maximum and the minimum 
densities for 50 protons. Following Ref.|l2| we denote the 
cells like a nucleus with Z protons and N neutrons. Thus, 
the two cells with Z=50 will be denoted by 1800 Sn and 
950 Sn. 

We consider also two representative cells with Z=40 
protons, namely 1500 Zr and 500 Zr. These cells corre- 
sponds to the densities p=0.0475 fm -3 and ,0=0.00159 
fm -3 , respectively. The cell 1500 Zr is the highest density 
cell with Z=40 while the cell 500 Zr corresponds to the 
lowest positive binding energy per nucleon. 

For the four cells chosen above we present first the pre- 
dictions of the HFB calculations for the densities and the 
mean fields of the nuclear clusters and the surrounding 
neutron gas. Then, we analyse how the distorsions of the 
density distributions induced by the clusters affect the 
pairing properties of nuclear matter inside the cells. 



A. Density and mean field distributions 

In finite nuclei there is always a maximum number of 
neutrons which can be bound for a given number of pro- 
tons. This neutron stability limit, which defines the neu- 
tron dripline, has been extensively studied theoretically 
in the recent years. Experimentally there are drastic 
limitations for approaching the neutron dripline in the 
laboratory since the neutron-rich nuclei are quickly beta 
decaying. This is not the case for the neutron-rich nuclei 
immersed in the inner crust since here the beta decay is 
blocked by the presence of the degenerate electron gas 
uniformly distributed throughout the baryonic matter. 
Consequently, inside the inner crust the nuclei can bind 
more neutrons than the nuclei in the vacuum. In addi- 
tion, their density and mean field can change significantly 
due to the presence of the surrounding neutron gas. 

These changes will be analysed first for the case of the 
Wigner-Seitz cells having 50 protons. The HFB calcu- 



lations with the SLy4 force and the pairing interaction 
mentioned in Section II predict for the isolated Sn iso- 
topes the dripline for two-neutron separation at N=126. 
Adding two more neutrons to 176 Sn the nuclear system 
in the vacuum is losing about 1 MeV from its binding. 




1800 Sn (n) 

— 950 Sn (n) 

176 Sn (n) 

--- 1800 Sn (p) 
---- 176 Sn (p) 



5 10 15 20 

r[fm] 

FIG. 1: Neutron and proton densities for the Wigner-Seitz 
cells with 50 protons and for the dripline nucleus 176 Sn. The 
cells to which the neutron densities (n) and proton densities 
(p) correspond are denoted by the total number of nucleons 
in the cell. 

In Figs. 1-2 the densities and mean fields of the dripline 
nucleus 1 Sn are compared to the corresponding quan- 
tities of the nuclear cluster immersed in the neutron gas. 
For the cell 950 Sn one can see that the nuclear cluster and 
the dripline nucleus 176 Sn have similar properties. More- 
over, the number of the bound neutrons inside the cluster 
is also equal to N=126. The situation is rather different 
for the cell 1800 Sn where the density of the outer neutron 
gas is much higher. Thus, it can be seen that the neu- 
tron density and the neutron mean field profiles have an 
extended " surface" before they reach the constant values 
corresponding to the neutron gas. In spite of these sig- 
nificant modifications of the density and the mean field, 
the number of bound neutrons inside the cluster, i.e., the 
neutrons with single-particle energies below the constant 
value of the mean field corresponding to the neutron gas, 
is still equal to N=126. On the other hand, if one inte- 
grates the neutron density up to r = 12 fm, where ap- 
proximatively the density profile is becoming constant, 
one finds about 300 neutrons. Thus, it appears that the 
large surface is generated in fact by the unbound neu- 
trons which are partially localised in the cluster region 
due to the interaction with the bound nucleons. 

Next, we discuss the properties of the nuclear clusters 
with 40 protons. The neutron dripline for isolated Zr 
isotopes is located at N=82. However, in this case the 
two-neutron separation energy curve is crossing zero with 
a very small slope. Compared to 122 Zr the binding en- 
ergy of 144 Zr is only about 200 keV smaller. This is due 
to the proximity of the resonant states 3/7/2 and 3^3/2 
to the continuum threshold and the fact that these states 
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FIG. 2: The mean fields for the Wigner-Seitz cells with 50 
protons and for the dripline nucleus 176 Sn. The notations are 
the same as in Figure 1. 



can easily accomodate more neutrons without losing too 
much binding energy [2l|. Thus, one expects that by 
immersing the dripline nucleus 122 Zr in the neutron gas 
one can easily gain the small extra energy necessary to 
bind more neutrons in the cluster. This is actually the 
case even for the low-density cell 500 Zr. The density and 
the mean field profiles for this cell arc shown in Figs. 3- 
4. From Fig. 3 one can see that although the density of 
the outer neutrons is about half compared to the cell 
900 Sn, they change more significantly the density pro- 
file of the dripline nucleus 122 Zr than for 176 Sn. Due to 
these changes the resonant states 3/7/2 and 3p 3 / 2 become 
bound, increasing by 12 the number of bound neutrons 
inside the nuclear cluster. 
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FIG. 3: Neutron and proton densities for the Wigner-Seitz 
cells with 40 protons and for the dripline nucleus 122 Zr. The 
notations are the same as in Figure 1. 



By increasing the density of the outer neutron gas one 
expects further changes of the nuclear cluster formed 



around 40 protons. These changes can be seen in Figs. 2- 
3 for the cell 1500 Zr. However, although the density of 
the outer neutron gas is now more than 20 times greater 
than in the cell 500 Zr, the number of bound neutrons in 
the cluster is increasing only by about 10 neutrons. The 
total number of neutrons found by integrating the den- 
sity up to r=12 fm is about 300 neutrons, as in the cell 
1800 Sn. 
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FIG. 4: The mean fields for the Wigner-Seitz cells with 40 
protons and for the dripline nucleus 122 Zr. The notations are 
the same as in Figure 1. 



In conclusion we find that the nuclear clusters im- 
mersed in the neutron gas keep many features of the 
dripline nucleus having the same number of protons. For 
the cells with large neutron densities one finds a thick 
surface-like structure which is developing in the transi- 
tion region from the center of the cell to the neutron gas 
regime. This surface structure is mainly formed by un- 
bound neutrons which are partially kept localised in the 
surface region by the interaction with the nucleons bound 
inside the nuclear cluster. 



B. Superfluid properties 

As we have seen in the previous subsection, the den- 
sity of the neutron matter is changing significantly from 
the center of the cell, where the nuclear cluster is lo- 
cated, to its edge filled by the uniform neutron gas. As 
shown by all the microscopic calculations, the pairing gap 
of the neutron matter depends strongly on the density. 
The majority of the calculations predict that as function 
of density the pairing gap has a bell shape form, with 
the maximum at a density equal to about one fifth the 
saturation density However, what is the maximum 
value of the gap and its detailed density dependence are 
still subjects of debate. Thus, the BCS-type calculations 
based on bare forces give a maximum gap of about 3 MeV 
and no pairing for nuclear matter at saturation density. 
Almost the same maximum gap one gets with a Gogny 
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force, but the gap falls to a finite value of about 1 MeV at 
the saturation density. On the other hand, all the micro- 
scopic calculations which go beyond the BCS approxima- 
tion show that the gap is strongly reduced by the screen- 
ing and self-energy effects. However, the amount of the 
gap supression and even the density dependence of the 
gap depend on the employed approximations. The most 
recent calculations point to a maximum gap of about 1 
MeV 13, llfl, |22f . How this value of the pairing gap could 
be reconciled with the pairing gap in finite nuclei is still 
an open question. 

Since at present it is not yet established what are the 
pairing properties of neutron matter, we perform here 
two calculations for the inner crust matter. In one calcu- 
lation we fix the density dependent delta force such as to 
reproduce approximately the gap values of neutron mat- 
ter provided by the Gogny force. The parameters of this 
pairing force, which we have used above for analysing the 
density and the mean field properties of the inner crust 
matter, are given at the end of Section II. In the second 
calculation, we fix the parameters of the pairing force so 
as to obtain for neutron matter a maximum gap of about 
1 MeV, as in the microscopic calculations which take into 
account the screening and the self-energy effects 0, ^3 ■ 
Since in this case the shape of the density dependence is 
unclear, we preferred to change only the strength of the 
force, which was reduced to the value Vo=-330 MeV fm 3 , 
and to keep the other two parameters as in the previ- 
ous calculations. For the cut-off energy we used in both 
calculations the same value, i.e., E c =60 MeV. 

First, we discuss the results given by the pairing in- 
teraction fitted to the Gogny force. The pairing field 
calculated for the two cells with 50 protons are shown in 
Fig. 5. For the cell 1800 Sn one can see that the pairing 
field is about two times smaller (here and below we refer 
to the absolute values of the pairing field) in the center 
of the cell compared to the value in the region of the uni- 
form neutron gas. We can also see that the pairing field 
is decreasing continuously towards the center of the cell. 

This is not the case for the cell 950 Sn. For this cell 
we observe that in passing from the low density region 
of the neutron gas towards the high density region of the 
cluster, the pairing field is increasing in the intermediate 
density region of the cluster surface. This is a manifesta- 
tion of the bell shape dependence of the pairing gap on 
density. Thus, for the cell 950 Sn the density of the outer 
neutron gas is much smaller than the value for which the 
gap is maximum, which is reached in the surface region. 
In the case of the cell 1800 Sn the density for which the 
gap is maximum is already reached in the neutron gas 
region. 

The behaviour of the pairing field in the cells with 
40 protons is rather similar with that found in the two 
cells discussed above. As seen in Fig. 6, for the cell 500 Zr 
we observe also an increase of the gap in the surface re- 
gion, which is now more pronounced than for the cell 
950 Sn. One can also notice that this significant increase of 
the pairing field in the surface region is extending rather 
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FIG. 5: The neutron pairing fields for the Wigner-Seitz cells 
with 50 protons. The full and the long-dashed lines corre- 
spond to the pairing fields calculated with the first pairing 
force while the dashed-dotted and the dashed lines correspond 
to the second pairing force. 



deeply towards the center of the cell. Consequently, the 
pairing field inside the nuclear cluster region is not be- 
coming anymore smaller than in the neutron gas region. 




1500 Zr(1) 

500 Zr(1) 

1500 Zr (2) 

500 Zr(2) 



12 



16 



20 



r[fm] 



FIG. 6: Same as Fig. 5, for Z—A0 protons. 

In all the cells we observe that the slope of the pair- 
ing field is changing very slowly when it is crossing the 
region between the nuclear cluster and the uniform neu- 
tron gas. This behaviour, sometimes referred to as the 
proximity effect Q. is related to the large size of Cooper 
pairs in low-density nuclear matter. The large diffusivity 
of the pairing field in the surface region, which is also 
seen in finite nuclei, is accentuated here by the fact that 
the superfluidity in the surface is carried mainly by the 
unbound neutrons, which are only partially localised at 
the interface between the nuclear cluster and the neutron 
gas. 

Next, we discuss the results corresponding to the sec- 
ond pairing force which simulates the polarisation ef- 
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fects. The pairing fields obtained with this pairing force 
are shown in Figs. 5-6 by short-dashed and dashed-dotted 
lines. As expected, the pairing correlations are reduced 
strongly for all the cells. For the high density cells 1800 Sn 
and 1800 Zr we can see that, in spite of the drastic supres- 
sion of the gap values, the profiles of the pairing fields are 
not changing very much compared to the previous calcu- 
lations. For the low-density cell 950 Sn we can see that the 
pairing field is becoming very small almost everywhere in 
the cell. However, in the cell 500 Zr the pairing field is still 
showing a significant increase in the cluster surface and 
it is almost vanishing in the neutron gas region. 

Thus, the behaviour of the pairing field in the inner 
crust matter is rather complex. As we have seen above, 
the nuclear clusters can not only suppress but also en- 
hance the pairing field inside the inner crust matter. On 
the other hand, the magnitude of the pairing field inside 
the inner crust depends very strongly on the scenario 
used for the pairing properties of infinite neutron mat- 
ter. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have analysed the properties of the 
inner crust matter in a self-consistent HFB approach. 
In the calculations we used for the particle-hole force a 
Skyrme-type interaction, i.e., SLy4, which is able to de- 
scribe properly the basic properties of finite nuclei and 
infinite neutron matter. As pairing interaction we have 
chosen a density-dependent contact force. The param- 
eters of the force were fixed so as to obtain for infinite 



neutron matter the pairing properties predicted by two 
calculations, i.e., a BCS calculation based on Gogny force 
and a microscopic calculation based on induced interac- 
tions. 

The HFB calculations were performed for a set of 
Wigner-Seitz cells representative of the structure of the 
inner crust matter. First we have studied the density and 
the mean field distributions inside the cells and we have 
analysed the properties of the nuclear clusters formed in 
the center of the cells. Then, we have studied how the 
pairing field is modified inside the cells by the presence 
of the nuclear clusters. We have thus found that, for the 
cells corresponding to high baryonic densities the pairing 
field is generally supressed in the region of the nuclear 
cluster. This is not the case for the low-density cells 
where the pairing field is increasing significantly in the 
surface region of the cluster compared to the region of 
the uniform neutron gas. 

The behaviour of the pairing field inside the cells is 
rather similar for both pairing forces. However, the val- 
ues of the pairing gaps are supressed dramatically if the 
pairing force which simulates the polarisation effects is 
used in the HFB calculations. 

The complex behaviour of the pairing field discussed 
here have important consequences on the thermal prop- 
erties of the inner crust matter. These aspects will be 
analysed in a forthcoming paper |'2.'"{| . 
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